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Abstract 

Several definitions for the average local value and local variance of a quantum observable are examined 
and compared with their classical counterparts. An explicit way to construct an infinite number 
of these quantities is provided. It is found that different classical conditions may be satisfied by 
different definitions, but none of the quantum definitions examined is entirely consistent with all 
classical requirements. 
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Many simple concepts of standard classical statistical mechanics are not easy to translate into 
quantum mechanics. The average local value and the local spread of a dynamical variable belong 
to this group. For an ensemble of particles in one dimension characterized by a joint distribution 
F(q,p) of position q and momentum p the average local value of the dynamical variable a(q,p) is 
given by 

a(Q) = p"j^y / / a(q',p)5(q' - q)F(q',p) dq'dp, (1) 
where P(q) = / F(q,p) dp, and the local spread, or local variance, by 

<4 = ^(?) - [%)F , (2) 

where a 2 is the local average of the square of a, or second order "local moment". Of course an 
arbitrary power of a, b(q,p) = a n (q,p), is also a function of (q,p), so the expression (JJ), mutatis 
mutandis, gives also the local moments of arbitrary order, 

%) = 1^ ffb(q',P)6(q'-q)F{q',p)dq'dp (3) 



P(q) 

= ^(q) = p^y // * n ti,PW - Q)FW,P) dq'dp. (4) 

This observation may appear trivial at this point but it will become important when discussing the 
quantum case. Also, the following marginal, joint, and conditional probabilities can be defined, 

P{a) = JJs[a-a(q,p)]F(q,p)dqdp (5) 
P(a,q) = J J 5[a- a(q' ,p)]5(q - q)F(q' ,p)dq dp (6) 
P(a\q) = P(a,q)/P(q) (7) 

where the last equation is simply Bayes' rule. In terms of these later quantities a and cr^. take the 
form 

®{q) = J a P(a\q) da = —— J a P(a,q) da (8) 



and 

•2 



= — J (a- a) 2 1 '(a. q) da . (9) 



P(q) 

The total (or global) average is simply the "g-average" [i.e. an average over q of a g-dependent 
function weighted by P(q)] of the local averages, 

(a) = J J a(q,p)F(q,p)dqdp = J a(q)P(q)dq, (10) 



whereas the total variance takes the form, 



/ J (a — (a)) 2 P(a,q) dadq (11) 
= J a 2 alq P(q) dq + J (a(q) - (a)) 2 P(q) dq . (12) 

a 2 is the (/-average of the local variances plus the "g-variance" of the local averages. This is an 
appealing decomposition, and it seems reasonable to seek for definitions of quantum local moments 
that preserve its structure. 

In a recent publication in this journal, L. Cohen ]T[ has argued that it is natural to interpret the 
following quantities as the local value of the observable associated with the hermitian operator A 
and its local spread, 

■ tm 
^ - mi • 

The superscripts S and C will distinguish these quantities from others to be defined later (the reason 
for using different letters will become clear soon), and the subscripts TZ and X indicate "real" and 
"imaginary" parts. (This interpretation of (|13D has also been proposed in |MJ].) These quantities 



obey a relation with the form of Eq. (12|), 



< = I * 2 A\MQ)\ 2 dq+ (A s (q) - (A)) 2 mq)\ 2 dq, (15) 



and the analogy with the classical equation has been invoked to understand the significance of the 
two terms in fll5|) JE]. In this letter we shall explore how far the classical-quantum analogy goes. We 
shall see in particular that the definitions ( |13"D and (0) do not always lead to results consistent with 
the form of equations (0) and (H). To this end it is useful to introduce an operator that symmetrizes 
A and 5(q - q) = \q)(q\, 



A q = ±[A6(q-q)+S(q-q)A), (16) 



and express ( [13]) in terms of it as, 



A S (q) = ^. (17) 
QKQ) 

where g(q) = \ip(q) | 2 is the probability density and (A q ) = {ip\ A q \ip). Note that the expectation value 
of A q is a "local density of A" , i.e., a quantity that integrated over q provides the global average, 



j(^\A q \^) dq=(A). 



(We shall see below that other definitions also satisfy this condition.) For consistency, the definition 
for the local average should apply to any operator, and in particular to A 2 . The local average of A 2 
is accordingly given by 

^ S(q) = ^sr = mA2 ^ - q)+ 5( $ - q)A2] w • (19) 

Following (U), a quantum local variance for A is "naturally" defined as 

[°% 9 ] 8 = A^ S (q)-[A s (q)] 2 (20) 
.^\[A 2 5{q-q) + 5{q-q)A 2 ]\i,) - [A s {q)] 2 (21) 



where the superscript S reminds that only symmetrized operators are used. This definition for the 



local variance also satisfies the decomposition fll2]) , but it is different from (0). To see this in more 
detail let us write (0) as 

2 / , , ^, ,, \ 2 



W 2 



c ( {q\A\1>) \ 

Alqi I m ) 



' (i>\A8(q-q)A\i>)-[A S (q)} 



{ <<#> 



(22) 
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Clearly the two procedures imply two different interpretations of the local density for the square of 
A. In general, 

(ij\A6(q- q)A\^} ? l -^\[A 2 5(q- q) + 5(q- q)A 2 }^) . (23) 

However their integrals over q are both equal to (A 2 ). Eq. (|20|) provides a local variance consistent 
with the definition given for the local average and it is in agreement with the classical expressions 
( [12| ) and (0), but it is not semidefinite positive. Its literal interpretation as a variance is therefore 
impossible. This resembles the status of the Wigner function and other phase space quasi-distribution 
functions. While not interpretable as probability distributions they can be used to correctly evaluate 
expectation values and to investigate the classical limit in a classical-like phase space language. 

In fact the decompositions of the total variance compatible with ( |12|) and (|]) are infinite. They 
can be found by means of the phase space formalisms described within the general framework 
provided also by Cohen Each of these formalisms is associated with a particular function 

f(6,r) of auxiliary variables 6 and r. The density operator and the operator A are related, re- 
spectively, to a quasi distribution function F(q,p; [/]) and a phase space representation A(q,p; [/]) 
[both depend functionally on /; Note that A(q,p; [/]) is not necessarily equal to the classical func- 
tion a(q,p)} in such a way that the expectation value of A is given by the phase space integral, 



(A) = JJ A(q,p; [f])F(q,p; [/]) dqdp. These formalisms are also closely related to the "quantization 
rules" , that define mappings from phase space functions to operators. The rules are generally used 
to associate a quantum operator with the classical function a(q,p). 

Local values can be defined using any of these phase space formalisms, see the classical expression 
d), as 



A(q; [/]) = J A(q,p; [f])F(q,p; [/]) dp . (24) 

Their g-average is the global average, (A) = j A gdq, and the variance can be decomposed in agree- 
ment with (O) by writing 

A 2 (q,p; [f])F(q,p; [/]) 



dp 



g(q)dq- / A\q-[f\) g(q) dq 



I A 2 (q;[f})g(q)dq-(A)\ (25) 



where A 2 (q,p; [/]) is the phase space representation of A 2 . 

The first two terms in (|25| ) may be regarded as the g-average of the "local variance" 

a% q [f]^(q;[f])-A\q;[f]), (26) 

while the last two terms take the form of a g-variance of local averages. An important point is 
that A 2 (q,p; [/]) is in general different from the square A 2 (q,p; [/]) so it is not necessarily a positive 
function. Unless / is suitably chosen F is not positive either, so that, contrary to the classical local 
variance, and in spite of the "square" used in the notation, a\, q [/] is not necessarily positive. Of 
course for a given kernel /, cr^J/] can be positive for particular observables and, similarly, for a 
given observable a family of kernels will give a positive local variances. In the case of momentum, 
A = p, Cohen has noted that (|2BJ) is positive for a family of kernel functions / that lead precisely to 
the choices (|X3j) and ([14]) [0. In general, however, each observable requires a different analysis. 

We shall next elaborate on the momentum observable and its powers using two of these formalisms 
because of their close relation to the previously discussed choices and to classical expressions. The 
Rivier-Margeneau-Hill (RMH) formalism is obtained by taking / = cos (#7-/1/2) H§,[5],||. We shall 
denote the corresponding representation of the state, namely the Margeneau-Hill function, as F MH . 
The subscript MH will be also used for the local averages and variances, ([24]) and (p6|), defined 
within this framework. The corresponding mapping from phase space representatives to operators 
is the "Rivier rule" ||. When applied to a phase space function with the factorized form g(q)h(p) 
it gives the symmetrized operator [g(q)h(p) + h(p)g(q)]/2. Note that for an arbitrary product of 



phase space functions AB, this symmetrization is not always equal to the simple symmetrization rule 
(AB + BA)/2, where A and B are quantum operators assigned to A and B by some prescription (that 



could be in fact Rivier's rule), see e.g. JlO|. However, if A(q) and B(p) are, respectively, functions of q 



and p only, and the associated operators are A(q) and B(p), then the two symmetrization procedures 
agree. This implies in particular that p i%S {q) = ~pF Mh \q) and [Cp^] 5 = { a p\^\ MH ■ 

It has been argued that the interpretation of (|T3|) as the local average leads to the Margeneau- 
Hill function [f|]. This connection can indeed be made, but, remarkably, the local values obtained by 
means of the RMH formalism in general differ from fllTf). To understand how this may happen let us 
briefly review the derivation of the Margeneau-Hill function in Q . The idea is to define a conditional 
probability density P (p\q) making use of the characteristic function concept, i.e., by inverting 

G(t, q) = e^ S {q) = J e lTp P s {p\q) dp (27) 



P b (p\q) = — / e- lTp e"P (q) dr (2* 
2n 



Expanding the exponential in e tTp , there results a series where the coefficients are the local averages 
of the powers of p, p n . Using the coordinate representation for the momentum operator the Taylor 
series for ip(q ± hr) can be recognized, 

5 (ir) n s ip{q + hr) iP*{q-hT) 



^ (q) = E (?) = y ; ( , ' + ; ,* ( , ' ■ ( 2 9) 

n n\ 2ip(q) 2ip*(q) 

If in addition one defines a quasi-joint probability distribution by the product g{q)P s {p\q), i.e., 
by formally adopting the structure of Bayes' theorem, this quasi-probability turns out to be the 
Margeneau-Hill function, 

F MH (q,p) = {(p\iP)(iP\q)(q\p)]n. (30) 

But, as discussed before, in general A s ^ A MH . The symmetrized operator {A 2 ) q is equal to the 
one provided by Rivier's rule if A is only a function of p, but will differ in general. The use of the 
structure of Bayes' theorem in the derivation of (|30"D is the key point that explains this seeming 
contradiction. This is a theorem valid for actual probabilities, but P s {p\q) and F MH (q,p) are not. 
In fact they can be negative. 

Finally, the comparison with the classical equations satisfied by local average and the variance is 
completed here by studying the time dependence, in particular the equations for time derivatives of 
the local averages of the powers of p, i.e., the equations of hydrodynamics. In this context the Weyl- 
Wigner (WW) phase space formalism is privileged. This formalism is associated with the simplest 

6 



choice, / = 1 |TT|,|5|.|(| . The phase space representative of the state is the Wigner function, F w (q,p), 
and the corresponding rule is the Weyl quantization rule. A superscript W will denote the local 
quantities (24]) and ( ^6l) calculated with this formalism. 



Cohen's prescription, Eq. fll3"D , and the phase space formalisms RMH and WW lead to the same 
equation for the local density, namely the continuity equation, 

T < 31 » 

where p = p s = p MH = p w . For the first local moment, the Wigner function leads to an "equation 
of motion" with exactly the same form as the classical one |T^| , 

dp w pWdp w d y {q) i d{ Q {g)[alX) 

dt m dq dq mg(q) dq 

m being the mass, and V(q) the potential function. (Even though the form is equal, the numerical 



values differ in general between the classical and the quantum cases, and the local variance [<r : 



2 ]W 
p\q\ 



is not semi definite positive ||12|| .) However, in the other approaches the local variances are different 

1 

',M l ""* r " + If(q) 



°IX = KX H + ^-<^/ - '/)/' -i»Hr, -<,)-»(<! -<,)!>) (33) 

1 



= K\ q }° -j^-j(2pS(q-q)p-p8(q-q)-8(q-q)p) , (34) 

and there appear extra terms which are not present in the classical equation (for the RMH formalism 
they were studied by Sonego fll3"| ). At least for several model potentials studied in ref. |I4| the WW 



formalism is not only closer to classical mechanics formally in this context; it is also numerically closer. 
The local kinetic energy densities g(q)p 2W (q) / (2m) , g(q)p 2MH (q)/(2m), and (p5(q — q)p)/(2m) were 
evaluated, and the one provided by the WW formalism was clearly the closest (numerically) to the 
classical values |14j|. Determining the extent of this agreement is an interesting open question for a 



separate study. 

In summary, we have examined different definitions of "quantum local averages and variances", 
and their similarities and differences with various classical expressions. None of them satisfies all 
the relations valid classically. It is useful to maximize the agreement with classical mechanics for 
examining the classical limit, and to gain physical insight in certain applications, but it should 
be noted that different classical criteria are satisfied by different quantum definitions. For specific 
applications one of the definitions may turn out to be the most convenient, but in fact each of 
them contains a piece of information which is only partially analogous to the corresponding classical 
quantity. 
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